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ABSTRACT 

We use photometric observations of solar-type stars, made by the NASA 
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Kepler Mission, to conduct a statistical study of the impact of stellar surface 
activity on the detectability of solar-like oscillations. We find that the number 
of stars with detected oscillations fall significantly with increasing levels of ac- 
tivity. The results present strong evidence for the impact of magnetic activity 
on the properties of near-surface convection in the stars, which appears to in- 
hibit the amplitudes of the stochastically excited, intrinsically damped solar-like 
oscillations. 

Subject headings: stars: oscillations — stars: interiors — stars: late-type - stars: 
activity - stars: magnetic field 

1. Introduction 

Solar-type stars show "solar-like" acoustic oscillations that are intrinsically damped 
and stochastically excited by near-surface convection (e.g., Houdek et al. 1999; Christensen- 
Dalsgaard 2004; Samadi et al. 2007). It is now well established that magnetic structures 
in the solar photosphere are strong absorbers of acoustic (or p-mode) oscillations (Braun et 
al. 1987, 1988; Braun & Birch 2008). A strong magnetic field can diminish the turbulent 
velocities in a convectively unstable layer (e.g., Proctor & Weiss 1982; Cattaneo et al. 2003) 
and this can affect the driving of acoustic modes (e.g., Jacoutot et al. 2008). In solar- type 
stars, the presence of a fibral magnetic field (Gough & Thompson 1988; Goldreich et al. 
1991; Houdek et al. 2001) may become sufficiently strong to affect not only the properties of 
the p-mode propagation, but also the turbulence of the convection by reducing its magnitude 
with increasing stellar activity, thereby reducing the amplitudes of the oscillations. 

The amplitudes (i.e., the square root of the total powers) of solar p modes are observed 
to decrease with increasing levels of solar activity (Chaplin et al. 2000; Komm et al. 2000). 
The decrease observed from solar minimum to solar maximum is about 12.5 % for modes of 
low spherical degree, / (Chaplin et al. 2000; Gelly et al. 2002; Jimenez- Reyes et al. 2003), 
which are the modes that are detectable in observations of solar-type stars. Garcia et al. 
(2010) recently uncovered the first evidence for changes in p-mode amplitudes associated with 
a stellar activity cycle in another solar-type star, from CoRoT satellite data on HD49933. 

HD49933, and the other F-type stars observed for asteroseismology by CoRoT, have 
activity levels that are not disimilar to those of the active Sun, as discerned from levels 
of variability in the lightcurves arising from rotational modulation of starspots and active 
regions (Mosser et al. 2005; 2009a). The same is true for the F5 star Procyon, as measured 
in both radial velocity (Arentoft et al. 2008) and photometry (Huber et al., 2010). However, 
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the G-type dwarf HD175726, which was also observed by CoRoT (Mosser et al. 2009b), 
shows much higher levels of activity but barely detectable solar-like oscillations. Mosser et 
al. speculated that the lower-than-expected amplitudes might have resulted from suppression 
by high levels of intrinsic magnetic activity. Dall et al. (2010) made a simlar suggestion for 
the active G8 star EK Eri. Here, we use the unprecedented large ensemble of oscillating 
solar-type stars observed by the NASA Kepler Mission to search for evidence of this effect 
in a large number of stars. 

In addition to searching for exoplanets (Borucki et al. 2010; Koch et al. 2010), Kepler is 
providing large quantities of high-quality data for the asteroseismic investigation of stars, as 
part of the Kepler Asteroseismology Investigation (Gilliland et al. 2010a). Photometry of a 
subset of these stars is being made at a cadence that is rapid enough to allow investigations 
of oscillations in solar-type stars, where dominant periods are of the order of several minutes 
(Chaplin et al. 2010,2011a; Christensen-Dalsgaard et al. 2010; Metcalfe et al. 2010). During 
the first seven months of science operations just over 2000 stars were observed for one month 
each as part of an asteroseismic survey of the solar-type part of the color- magnitude diagram. 
Solar-like oscillations have been detected in about 500 stars, increasing by a factor of about 
25 the number of solar-type stars with detected oscillations. 

The large number of solar-type stars in this Kepler ensemble makes possible the statis- 
tical study of intrinsic stellar properties and trends, in what is a homogenous data sample of 
unprecedented quality. Here, we use results on the ensemble to conduct a statistical study 
of the impact of stellar activity on the detectability of solar-like oscillations. 

For this study we have used two simple measures of variability in the Kepler lightcurves 
as proxies of the levels of stellar magnetic activity, including one measure suggested by 
Basri et al. (2010, 2011) from their survey of activity levels in more than 100,000 stars 
observed in 30-min (long) cadence by Kepler during its first month of science operations. 
We study a subset of about 2000 solar-type stars having short-cadence Kepler data and test 
if the distribution of observed variability differs for stars that do, and do not, have detected 
solar-like oscillations. 



2. Data and Analysis 

We use asteroseismic results on solar-type stars that were observed by Kepler during 
the first seven months of science operations. About 2000 stars, down to Kepler apparent 
magnitude Kp ~ 12.5, were selected as potential solar-type targets based upon parameters 
in the Kepler Input Catalog (KIC; Batalha et al. 2010, Koch et al. 2010). Each star was 
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observed for one month at a time in short- cadence mode (58.85s; see Gilliland et al. 2010b). 
Time series were prepared for asteroseismic analysis in the manner described by Garcia et 
al. (2011), using procedures that work on the raw lightcurves. Lightcurves prepared for 
the Transiting Planet Search pipeline are not appropriate for use here since phenomena 
associated with stellar variability are suppressed or removed to aid planet detection (Jenkins 
et al. 2010). 

Different teams analyzed the prepared lightcurves to attempt to detect signatures of 
solar-like oscillations (see Huber et al. 2009; Mosser & Appourchaux 2009; Roxburgh 2009; 
Campante et al. 2010; Chaplin et al. 2010, 2011b; Hekker et al. 2010; Karoff et al. 
2010; Mathur et al. 2010). Most of the applied detection techniques relied on extracting 
signatures of the near-regular frequency separations of the solar-like oscillation frequency 
spectrum, while others searched for signatures of the Gaussian-like power excess due to the 
oscillations. Verner et al. (2011) present a detailed comparison of the results returned by the 
different pipelines on the Kepler lightcurves, finding reassuring levels of agreement between 
the results. We demanded that at least two of the asteroseismic data analysis pipelines 
returned consistent results on a star for it to be flagged as a solid detection for use in this 
paper. A total of around 500 stars were flagged as having detected solar-like oscillations. The 
signal-to-noise ratio in the oscillations required for a detection was > 0.1, as defined by the 
ratio of the maximum power spectral density of the smoothed oscillation envelope relative to 
the estimated background at the frequency of maximum oscillations power. This threshold 
level allows an unambigous detection of the signature of the large frequency separation, and 
identification of a significant power excess due to the oscillations. 

For each star that was analyzed, we used two simple metrics of variability as proxies of 
the stellar surface activity, as determined from direct analysis of the prepared lightcurves. 
For both, we first applied a low-pass filter to remove the oscillation signal smoothing each 
lightcurve with a one-hr-long boxcar. For one metric, we measured the maximum absolute 
deviation of each smoothed lightcurve from its mean. This simple measure of variability is 
what Basri et al. (2010) call the "range", and here we label it r^. Basri et al. smoothed 
long-cadence (29.4 min) Kepler lightcurves with a 10-hr boxcar. We tested the effect of 
smoothing our short-cadence lightcurves on timescales ranging from 1 hr up to 1 d, but 
found no significant impact on our results. As our second variability metric we use what 
Garcia et al. (2010) referred to as a "starspot proxy" . We measured the standard deviation 
about the mean (i.e., the RMS) of each smoothed lightcurve, and we call this metric <7\ ir . 

Neither of the measures, as defined above, takes explicit account of the apparent mag- 
nitude of the target. We would expect there to be a magnitude-dependent correction to r hr 
and <Thr on account of the changing contribution due to shot noise, although the effect is 
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small. We specified the required corrections using the "minimal noise" model for Kepler in 
Gilliland et al. (2010b). The RMS noise, a, per At = 58.85-sec integration is given by 

a = — (c + 9.5 x 10 5 (U/Kp) 5 ) 1/2 ppt, (1) 
c v 

where c = 1.28 x io°- 4 ( 12_i:s: P)+ 7 detections per cadence, and Kp is the Kepler apparent 
magnitude. Since the time series are smoothed with a 1-hr (3600-sec) boxcar, the additive 
correction that must be removed from o"h r is just a(At/3Q00) 1 ^ 2 . The correction for r^ T was 
calibrated with Monte-Carlo simulations, and found to be ~ (2/3)<r. 

Since our aim is to understand the impact of stellar activity on the detectability of the 
solar-like oscillations, we ignored those stars whose lightcurve variability could be attributed 
to another phenomenon, e.g., eclipsing binaries, and classical pulsators at the hot end of the 
sample. The fraction of eclipsing binaries that we removed (~ 1 %) is in line with the rate of 
occurence of binaries in the Kepler field of view, as determined by Prsa et al. (2011) down 
to Kp ~ 16. 



3. Results 

The top two panels of Fig. CD plot the range, r^ r , and the RMS, ov, as a function of 
T e g. Points in black are stars with detected solar-like oscillations. The bottom two panels 
plot the metrics with the Kepler apparent magnitude, Kp, as the independent variable. The 
dotted lines follow the additive corrections (see above). 

Our range plot in Fig. [T]is very similar to the corresponding plot in Basri et al. (2010), 
which shows results on just over 100,000 stars observed during the first full month of Kepler 
science operations. The data have a lower-limit envelope that has its minimum at approxi- 
mately solar temperature. The lower limit of the envelope shifts to higher levels of variability 
at lower and higher temperatures. 

It is also apparent from Fig. [T] that the cloud of points for stars with detected solar-like 
oscillations (black points) is clustered toward lower levels of variability. We see no detections 
above r^ r of ~ 20 ppt, and of ~ 10 ppt. The cloud of points for stars without detected 
oscillations (gray points) extends to much higher values in both r^ r and <7h r - There are also 
stars having no detections that nevertheless show the same levels of variability as stars with 
detections, i.e., at the lower levels of variability the distributions overlap. The apparent 
gap at T fj ~ 5300 K in the distribution of stars with detected oscillations is the result of 
higher numbers of detections at lower T e ff from evolved stars at the base of the red giant 
branch, which have higher oscillation amplitudes than their main-sequence cousins. The 
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Fig. 1. — Range, 7V, and RMS, o"hr as a function of T e g (top panels) and Kepler apparent 
magnitude, Kp (bottom panels). Stars with detected solar-like oscillations are plotted in 
black; stars with no detections are plotted in gray. The dotted lines follow the additive 
corrections that were applied to r hr and cr hT (see text). 
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Fig. 2. — Top panels: Histograms of r hr and <7 hr , for all analyzed stars (gray solid lines) and 
stars with detected solar-like oscillations (black dashed lines). Middle panels: Histograms 
normalized to a maximum value of unity (same linestyles). Bottom panels: Fraction of stars 
in each histogram bin showing detected solar-like oscillations. 
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distributions in Kp (bottom panels of Fig. [TJ) are also not smooth, but this effect is present 
for stars both with, and without, detected oscillations, and it is merely a selection effect 
arising from the final choice of target stars. 

Histograms of and a^ r are plotted in the top panels of Fig. El for all analyzed 
stars (gray solid lines) and stars with detected oscillations only (black dashed lines). The 
middle panels show the histograms normalized to a maximum value of unity to allow a visual 
comparison to be made of the shapes of the distributions of stars with, and without, detected 
oscillations. The bottom panels plot the fraction of stars in each histogram bin that show 
detections, with the errors calculated from Poisson statistics. The histograms in the middle 
panels show the aforementioned significant deficit of detections at higher levels of variability. 
The significance of this fall-off is confirmed by the detection ratios in the bottom panels. 

Two categories of explanation for the fall-off suggest themselves: one that has its origins 
in the intrinsic properties of the stars (much the more interesting explanation); or another 
that is the result of either data-analysis or selection bias issues affecting the detectability of 
the modes. Let us consider first the possibility that data analysis issues might be the cause. 

Sudden discontinuities or large excursions in the lightcurves can affect the appearance of 
the frequency spectrum of the lightcurves if not treated properly, e.g., as a result of the intro- 
duction of a complex overtone structure, or the introduction of frequency-dependent noise, 
into the frequency range of interest for detection of the oscillations. To test the impact of 
these effects on our results we selected about 150 stars for which we had successfully detected 
oscillations, covering a wide range of intrinsic stellar properties and apparent magnitudes. 
We then picked a subset of 50 stars with high levels of variability, but no detected oscilla- 
tions. In one set of tests - analyzed using the pipeline described by Hekker et al. (2010) - we 
added the low-pass filtered parts of the high variability lightcurves to the high-pass filtered 
parts of the lightcurves with detected oscillations, and then checked whether we could still 
detect the oscillations. In another set of tests - analyzed using the pipeline described by 
Huber et al. (2009) - we selected high variability stars at random, and generated synthetic 
lightcurves comprised of 10 random low-frequency sinusoids with amplitudes selected to en- 
sure that the low-frequency spectra of the synthetic lightcurves matched those of the real 
stars. These synthetic lightcurves were then multiplied by a random factor - to allow as 
wide a range of variability to be sampled as possible - and then added to one of the selected 
lightcurves with detected oscillations. Again, we tested to see whether the oscillations could 
still be detected. The data-analysis pipelines detected oscillations at all levels of variability 
(i.e., even in the most extreme cases), and while a modest fall-off of the detection rates was 
seen at the highest variability (e.g., for > 30ppt) this was at nothing like the levels seen 
in the real Kepler data. 
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Another concern would be that the fall-off is the result of selection bias: for example, 
that stars showing higher levels of variability also tend to be fainter, on the average, and 
therefore less likely to show detected oscillations due to higher levels of shot noise. There 
are certainly more stars at fainter Kp, which we would expect to more fully sample the 
underlying distribution of r^ r and o"h r , including the higher values. The bottom panels of 
Fig. [1] do seem to bear this out. To test the possible effects of selection bias, we ran the 
detection prediction code developed for use by the Kepler Science Team (Chaplin et al. 
2011b) on all 2000 targets. This code takes as input the KIC radius, effective temperature 
and apparent magnitude, and produces as output an estimate of the probability of detection 
for an assumed length of observation, based on use of scaling relations that know nothing 
about the possible effects of activity on the amplitudes of the oscillatons. We are interested 
only in the distribution of predicted detections, and the ensemble is large enough to provide 
statistically robust results. Those results indicated that the steep fall-off seen in the Kepler 
results cannot be explained by selection bias. There was some fall-off in the predicted fraction 
of detected oscillations, because there are more stars with high variability that are also faint 
(see above); but we saw "predicted" detections in stars showing even the highest levels of 
variability. The results also provided an explanation of why some stars with low levels of 
variability did not have detected oscillations. The expected success rate of detections was 
found to be less than 100 %, because noise realizations will in some cases hamper extraction 
of the oscillation signals. That said, the absence of detections in some of the brightest 
stars in the sample that show only modest levels of variability is a puzzle. It may be that 
the inclinations of stars play a role here. The inclination will affect the apparent (i.e., 
observed) variations seen in r hr and o"h r , assuming that in solar- type stars those variations 
are dominated by contributions from active latitudes like for the Sun (e.g., see Knaack et al. 
2001; Vazquez Ramfo et al. 2011). It could be that some stars with small tv and o"h r are 
actually intrinsically active stars observed at low angles of inclination (which will reduce the 
variations observed in the lightcurves). 

In summary, we conclude that the steep fall-off in the observed fraction of detections 
has a stellar explanation. 

4. Discussion 

The most compelling explanation for the results is that they show evidence for intrinsic 
stellar (magnetic) activity suppressing the amplitudes of the solar-like oscillations, and hence 
adversely affecting the detectability of those oscillations. In offering this explanation we 
assume that the metrics we have used are reasonable proxies of intrinsic levels of stellar 
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activity (see Basri et al. (2010) for further discussion). This certainly seems to be the case 
for the Sun. We analyzed observations of the bolometric flux of the Sun made by the PM06 
instrument onboard the ESA/NASA SOHO spacecraft (Frohlich et al. 1997). Fig. [3]plots r hr 
and ov (filled symbols in each panel) as determined by analysis of independent one-month- 
long segments of PM06 data. The gray line in each panel describes a smooth curve through 
the independent measures, while the dotted line shows the scaled 10.7-cm radio flux, which 
is an excellent proxy of the global level of magnetic activity on the Sun (e.g., see Chaplin et 
al. 2007, and references therein). The underlying trends in tv and a^ r clearly follow those 
in the magnetic activity. 

The range varies from about 0.25 ppt to about lppt between solar minimum and max- 
imum. We know that the amplitudes of the low-degree solar p modes are at the same time 
suppressed by a fraction ~ 0.125 (Chaplin et al. 2000; Gelly et al. 2002; Jimenez-Reyes et 
al. 2003; Garcia et al. 2010). These solar data in principle allow us to calibrate the expected 
suppression of oscillation amplitudes due to activity. Here we present a simple estimate for 

We assume that the fractional amount by which the amplitudes are suppressed, SA/A, 
is a linear function in r^. This is a valid assumption for the solar values, but may of course 
be questionable at higher levels of activity. There is also the impact of the inclination of the 
star to consider (see above). With these caveats in mind, we have that 



Integration of the above gives the resulting suppressed amplitude, expressed as a fraction of 
the amplitude expected for no activity (i.e., zero rv): 



Fig. Hlplots A(rhr)/^4(0) as a function of r^. At a range of 20ppt, the prediction is that the 
amplitudes are suppressed by a factor of almost 30, and it is therefore not surprising that 
we see hardly any detections in the Kepler ensemble at this value, and none above. 

The predictions in Fig. H] are also in agreement with the CoRoT results on the active 
G-type dwarf HD175726 reported by Mosser et al. (2009b). Peak-to-peak variations in 
the lightcurve, due to rotational modulation by spots, were found to be typically 1 %, so 
that TV = 5 ppt. Mosser et al. measured amplitudes that were about 1.7-times lower than 
expected, based on scaling- relation predictions that are calibrated against solar (i.e., low 
activity) values. Fig. H] implies that at this r^ T amplitudes should be suppressed by a factor 
of about two, very close to the factor reported by Mosser et al. 




(2) 



A( rhl )/A(0) ~ exp (-r hr /6) . 



(3) 
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Around 100 of the Kepler stars showing detected solar-like oscillations will be observed 
for periods lasting several months up to a few years. This will in principle allow us to further 
constrain the effects of magnetic activity on the oscillation amplitudes by measuring changes 
to the amplitudes, and the resulting detectability of the modes, as activity levels vary in 
time. An additional 100 stars, selected by the Kepler Science Team as possible planet hosts, 
should show oscillations based on Chaplin et al. (2011b) and will be followed for several 
months to years. 

Funding for this Discovery mission is provided by NASA's Science Mission Directorate. 
The authors wish to thank the entire Kepler team, without whom these results would not be 
possible. We also thank all funding councils and agencies that have supported the activities 
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Fig. 3. — Estimates of r hr (left-hand panel) and a^r (right-hand panel) for the Sun (filled 
symbols), as determined from analysis of one- month- long segments of PM06 data. Gray 
lines describe smooth curves through the independent measures, while the dotted lines show 
the scaled 10.7-cm radio flux. 
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Fig. 4. — A simple model of the suppression of mode amplitudes by stellar activity, with 
activity measured by the range parameter, r hr . The figure shows the expected amplitude 
versus r hr , as a fraction of the amplitude expected for zero activity (i.e., r hr = 0). 



